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This paper theoretically discusses the photon emission spectra of a superconducting pn-junction. 
On the basis of the second order perturbation theory for electron-photon interaction, we show 
that the recombination of a Cooper with two p-type carriers causes drastic enhancement of the 
luminescence intensity. The calculated results of photon emission spectra explain characteristic 
features of observed signal in an recent experiment. Our results indicate high functionalities of 
superconducting light-emitting devices. 

PACS numbers: 74.50.-|-r, 74.25.Fy,74.70.Tx 
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Light-emitting diode (LED) usually fabricated on 
semiconductors has been an important element of mod- 
ern technologies. A trend of research seems to be focusing 
on producing a better controlled photon and an entan- 
gled photon pair for realizing quantum computa- 
tion and quantum information. Superconducting devices 
have an advantage to obtain robustly coherent quantum 
states because of its coherent nature 0, 0, @ • Supercon- 
ducting LEDs [6] have been originally proposed in the 
context of superradiation. They are, however, a promis- 
ing candidate to create an entangled photon pair f?]. A 
recent theoretical study predicts the Josephson radiation 
in a superconducting pn junction Thus supercon- 
ductor/semiconductor LED hybrids undoubtedly have a 
possibility to produce technologies in the next generation. 

The radiative recombination of Cooper pairs has been 
observed recently in a InGaAs/InP pn junction attaching 
onto a superconductor Nb [8]. The electroluminescence 
becomes drastically large at low temperatures below the 
superconducting transition temperature Tc of Nb elec- 
trode. Surprisingly degree of the enhancement in the lu- 
minescence intensity is one order of magnitude. Although 
the effects of superconductivity on the radiative recom- 
bination are clear in experiments, a mechanism has been 
an open problem. We theoretically address this issue in 
this paper. We study the emission spectra of photon in a 
superconducting pn junction based on the second order 
perturbation theory. In the second order expansion, we 
find that a peculiar recombination process to supercon- 
ductivity enlarges the luminescence intensity, where two 
electrons recombine with two p-type carriers as a Cooper. 
The theoretical results explain characteristic features of 
the experimental findings [Sj . This paper not only figures 
out a mechanism of the large amplitude of luminescence 
intensity but also gives a guide for designing highly func- 
tional superconducting light-emitting devices. 

Let us consider a p-type semiconductor / supercon- 
ductor junction as shown in Fig. [TJa). The energy is 
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FIG. 1: (Color online) Schematic energy diagram of pn junc- 
tions. A theoretical model used for calculation is shown in 
(a). In (c), a real junction in an experiment is illustrated. 
Theoretical results of the photon spectra are shown in (b) 
and (d). 



measured from the horizontal line indicated by '0'. The 
sign of energy in a p-type semiconductor is chosen to be 
opposite to that in a superconductor. We assume that 
a semiconductor and a superconductor are in their local 
equilibrium which are characterized by the local chemical 
potential /ip and /i„ , respectively. The edges of the con- 
duction and valence bands are Ec and Ey, respectively. 
In what follows, we use a unit oi h — kB = c = \, where 
kB is the Boltzmann constant and c is the speed of light. 
The p-type semiconductor is described by 



fc,fT 
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where ep(fc) — A;^/(2mp)-f i?t, + el^sd/2, rup is the effective 
mass, Vsd is the applied bias voltage across the junction. 
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^^'^ a(^k a) Creation (annihilation) operator of 

a p-type carrier with a wave number k and spin a or 
|. The photon states is described by 



(2) 



where a^(ag) is the creation (annihilation) operator of a 
photon with a wave number q and an energy ujq. The 
normal state in a metal is described by 



k,< 



with m„ being the effective mass. The electron-photon 
interaction Hamiltonian reads 

-ff / = ^ Bk^q bk-q,cr Ck.a a\ + H.C., (4) 
k,q,a 

where Bk q is the coupling energy. On the basis of the 
second order perturbation theory, the number of photon 



N„h = J2q '^q'^q '^^ calculatcd as 



(5) 



(Nphil)) = dh dt2{xo\Hi{h)NphHi{t2)\xa), (6) 



dU 



dts / dU I (2), 
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I{2) = {xo\Hi{h)Hi{t2)NphHi{t:i)Hi{U)\Xo), (8) 
Ixo) -> |0)® 1^)® |P), (9) 



where |0) is the zero photon state. 

The BCS theory describes superconducting states, 



fe,(T 
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where Eu = V^IW+A^, ^n{k) = fcV2m„-/i„, A is the 
pair potential, and 7^0.(7^0-) creation (annihila- 

tion) operator of Bogoliubov quasiparticle. This descrip- 
tion, however, is valid within a small energy scale near 
the Fermi level which is at /i„ = i?c + eVsd/2 + mea- 
sured from '0'. To apply the BCS theory to the present 
issue, a rule is necessary to describe the operator in the 
interaction picture. The canonical transformation con- 
nects an electron operator and Bogoliubov operators by 



s.Vke-'^-'-f_^-\ , (11) 



in (xol • • • Ixo), where uM = [{1 + {-)^n{k) I E^) 12^/^ 
So- — 1(— 1) for a =t (i), and a means the opposite spin 
to a. The thermal average of operators is carried out 
in the local equilibrium. In a p-type semiconductor, for 
instance, the average of operators are calculated in 



(12) 



instead of Eq. ([T]) with ^p(fc) = k^/2nip — /Xp. In a super- 
conductor, the average of the Bogoliubov operators are 
calculated in Eq. (fTO|l . In Eq. ([9]), \P) means the state 
vector of p-type carrier in the local equilibrium and \N) 
indicates the BCS state in the local equilibrium. 

The time average of the photon number {Nph) corre- 
sponds to the luminescence intensity and it in the first 
order perturbation expansion results in 



(iVp,(l)) =27t}_^ \Bk.qVfLq 



X [ulfl!5iCu - Ek) + vl{l - r^)5{Cj + Ek)] , (13) 



where Co = 



Eg fin flp Cp(^ 



Eg = 



E,+E, + eVsd. fl! = [1 - tanh(i;fc/2T)]/2, and = 
[1 — tanh(^j,(fc)/2r)]/2. This result recovers the photon 
spectra in a normal pn junction by tuning A ^ 0, which 
means that Ek — > —inik),Uk ^ 0, Wfe — > 1 for k < kp 
and Ek £.nik),Uk ^ 1, Wfc ^ for fc > kp with kp 
being the Fermi wave number satisfying kp/2'mn — /i„. 
The threshold of spectra is Eg and the width of spec- 
tra is given by fin + fip- The spectra in Eq. (jl3p have a 
broad profile reflecting the quasiparticle density of states 
as shown in Fig. [Tfb). 

One of the characteristics peculiar to superconductiv- 
ity is a giant oscillator strength This comes from 
the much freedom of the wave vector k for the remaining 
elementary excitation when a Cooper pair decays radia- 
tively with a p-type carrier with the opposite wave vector 
~k. This would gives the large luminescence intensity 
in the experiment Q . This effect of the giant oscillator 
strength works also in the first stage in the second order 
expansion while not in the second process. The other 
characteristics is coming from the resonant enhancement 
due to the nearly degenerate nature of the final and in- 
termediate states close to the initial state. This is shown 
in the following calculation. 

The results of the second order perturbation are given 

by 



m = E 



fci---fc4,qi---q4,o"i---o"4 

X Bi^.q^Bl^ q^Bk^^q^Bk^^q^QphQpQN: (14) 

with Qj{kj, Qj) = LUq. — ep{kj — g^) — fj-n- The average of 
the operators Qph, Qp and Qn are calculated as follows. 



Qph =E<0| W"lv«|o>' 



la 

=2{5'i,5l^ + 51,51,), (15) 

where 51^ = 5q.^q., Jf^ = 5^^,^^, 5f^ = 5p^^p. and = 
kj — Qy By applying the Bogoliubov transformation, we 
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find, 



'''''^''iU..,.m. (17) 



which gives twelve terms. In what follows, we extract the 
most dominant contribution in the second order terms. 
The average of Q at includes following four terms 

Qn{S) — UkiVkiUk3Vk3Sai,a2Scr3,94^ki,-k2Sk3,-ki<^l<^3 

X gi^ki (tl-t2)g-i-Efc3(t3-t4) jn (1 — ) 

_ g-i-Bfci(tl-t2)g-iBfc3(t3-t4)^]^ _ j-n ^^j^ _ j-n 



(18) 



Because da-3,g-i5k^^^k^ in Eq. (|18p means the destruction 
of two electrons as a Cooper pair, Qn{S) describe effects 
of superconductivity on the emission spectra. Another 
eight terms in Qn contribute to the emitting processes 
described by Fig. [H[b) which gives the luminescence in- 
tensity proportional to {Nph{l))'^. We will show that 
Qn{S) gives a large contribution to the emission spec- 
tra at 6q = Qi + q2 = 0, {Sq = q^ + q^ = in other 
wards) Substituting Eqs. (US]) , (HH) and HI]) into 

Eq. p4p and carrying out time integrations, we obtain 



(19) 
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[Ek ~ ijrf {Ek + i/t) 
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(20) 



where we introduce a relaxation time t to remove effects 
of the perturbation at i — s- —00, we neglect dependence of 
B on wave numbers and assume f^^_q = 1. At l/r = 0, 
/o = 7rAfo/2A essentially diverges for small A with A^o 
denoting the normal density of states in a superconductor 
at the Fermi energy. The singular behavior at small A 
in Eq. (|20p is a sign of the large luminescence intensity 
due to superconductivity. 

We first show mathematical reasons of the singular- 
ity. Then we will discuss the physics behind the phe- 
nomenon. A two-photon emitting process in Qn{S) is 
illustrated in Fig. [Ha). The annihilation of a Cooper 
pair is described by c_fc^|Cfc | which includes a opera- 
tor 7^ ^7^ |. Let us assume that the energy of the ini- 
tial state is zero. In the first order expansion, the op- 
eration of 7fcj to the BSC state decreases energy by 
Ek + An- At the same time, a p-type carrier with en- 
ergy ep(fc — q) is destructed and a photon with energy Uq 
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FIG. 2: (Color online) Recombination processes in the second 
order perturbation expansion, where solid, broken and wavy 
lines represent the propagation of an electron, a p-type carrier 
and a photon, respectively In (a), a recombination of a Cooper 
pair in Qm{S) is shown. In (b), a recombination process other 
than Qn{S) is illustrated. 



is created. Thus the energy of the intermediate state 6E1 
results in SEi ^ ujq ~ ep{k - q) - Ek - fin = ^k,q - Ek 
which is the energy denominator in the perturbation ex- 
pansion. In the second order, the operation of 7^ |, 
the destruction of a p-type carrier, and the creation of 
a photon gain energy by Ek — fin, —^p{—k + q), and 
Lo^q, respectively. Therefore the difference in energy be- 
tween the intermediate state and the final one becomes 
5E2 = ujq—ep{k—q)+Ek—fLn = flk,q+Ek- The perturba- 
tion theory requires the energy conservation between the 
initial and the final states, (i.e., SEi + SE2 ~ 0), which 
leads to 2ilk q ~ 0. As a result, only a small value of Ek 
remains in the denominator as shown in Eq. (|20p . The 
physics behind the phenomena is simple. The BCS state 
can have ability to emit a pair of photons with remain- 
ing its state almost unchanged because the BCS state is 
the eigen state of 7^ ^7^ ^ The equation ^kp.q = is 
the condition for emitting a photon. The threshold and 
width of spectra are Eg + /i„ and fj,p, respectively. In 
Fig. [Ijb) , we show predicted spectra in the second order 
process. 

The singular behavior in perturbation expansion im- 
plies an importance of higher order terms to predict 
the luminescence intensity quantitatively. Here we do 
not discuss this issue, but choose an alternative way 
of regularizing the obtained results for qualitative argu- 
ment. In what follows, we introduce a finite relaxation 
time. First we consider mean free time due to elastic 
impurity scatterings tq. At T = 0, we obtain Iq = 
Ioo{0)2a'^/{\/TTa^{a + Vl + a^)), where a = tqAq, 
Ao is the pair potential at the zero temperature and 



/oo(0) = 7rA^o/2Ao is Eq. ([20]) at T 
we find 



0. At T < Tc 



coa^{A/AofAo/T a<l 
/oo(0) ' I a^(A/Ao)2 a» 1, 



(21) 
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where Cq is a constant of the order of unity. In Fig. [3lja), 
we show Iq as a function of temperature for several 
choices of a, where we describe the dependence of A on 
temperature by the BCS theory. The amphtude of Iq at 
T = is suppressed in the dirty hmit as shown in a re- 
sult with a = 0.2. The amplitude at T = increases with 
increasing a. At a = 1, /o(0) has almost the same ampli- 
tude as /oo(0). When we increase a up to 2.0, the results 
show a bump just below Tc. Next we consider inelas- 
tic scatterings described by l/ne — de [T/Tcf , where 
Cie is a coupling constant and p depends on scattering 
sources such as p = 1 for electron-phonon scatterings 
and p = 2 for repulsive electron-electron interaction. In 
Fig-IHb), we calculate Jo for several choices of Cie and p. 
Since l/xie — > at T = 0, the amplitude is close to /oo(0) 
at T = 0. When we decreases Cie, the bump appears be- 
low Tc as well as in (a). For l/r < Aq, the luminescence 
intensity at T < Tc is then given by 



{Nph{2)) ~ AttcoIBI^Nq 



T 



(22) 




FIG. 3; Temperature dependence of luminescence intensity. 
In (a), we consider relaxation time due to the elastic impurity 
scatterings hy a — tqAo. In (b), we consider the relaxation 
due to inelastic scatterings. 



Finally we modify Eq. ([22|) to describe the photon spec- 
tra in the experiment 8] as shown in Fig. [TJc). In the 
real junction, a superconductor is attached to a n-type 
semiconductor whose thickness is L^. The proximity ef- 
fect enhances the luminescence intensity. In Eq. (P^ . 
A is proportional to the amplitude of a Cooper pair. 
In n-type semiconductor, the proximity effect enables 
the pair amplitude which proportional to Ae~^™/^^ with 
= y/D/2TTT and D being the diffusion constant in the 
n-type semiconductor. The photon pairs arc emitted in a 
quantum well between the p- and n-type semiconductors. 
The level in the quantum well should coincide with 
the Fermi level in the n-type semiconductor Namely 
\Eii] — fin\ must be less than both the Thouless energy 
Exh — D/L"^ and A. This resonant condition is partic- 
ularly important for a Cooper pair to penetrate into the 



quantum well. The emission spectra has a peak at ujq 
and the peak width is given by F = t^No, where is 
the transfer integral between the quantum well and the 
semiconductor. The argument above is summarized by 
an equation for T <Tc 



{N,nm)^\B\^N,TY, 



K-o.oP + (F)2' 



(23) 



where we introduce the Lorentz resonant function by 
hand. In the experiment, is estimated to be much 
larger than below Tc. Thus the theoretical results 
in Fig. [3] can describe experimental results of the lumi- 
nescence intensity. In fact, the experimental results of 
Fig. 6(b) in Ref. |8| show a very similar line shape to that 
in Fig.[3Ka) with a — 1. 

In conclusion, we have studied the photon emission 
spectra in a superconducting pn junction based on the 
second order perturbation theory for electron-photon in- 
teraction. We have found in the second order expansion 
that a peculiar recombination process to superconductiv- 
ity enlarges the luminescence intensity. The theoretical 
results explain temperature dependence of the lumines- 
cence intensity observed in an recent experiment. 
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The speed of light is much larger than the Fermi velocity. 
In such case, this condition becomes |qj^| ~ Iqjl- 



